Theory of spin relaxation in magnetic resonance force microscopy 
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We study relaxation of a spin in magnetic resonance force microscopy (MRFM) experiments. We 
evaluate the relaxation rate for the spin caused by high-frequency mechanical noise of the cantilever 
under the conditions of adiabatic spin inversion. We find qualitative agreement between the obtained 
relaxation time and the experimental results of Stipe et al. [Phys. Rev. Lett. 87, 277602 (2001)]. 
Based on our analysis, we propose a method for improving the MRFM sensitivity by engineering 
cantilevers with reduced tip positional fluctuations. 



Magnetic resonance force microscopy (MRFM) has 
proved to be a powerful tool in studying magnetic prop- 
erties of materials [1-4]. Unlike conventional magnetic 
resonance techniques, MRFM allows one to probe magne- 
tization with a nanometer-scale spatial resolution, which 
is important for practical realization of spintronic [5] and 
quantum information processing devices [6]. Measure- 
ment of magnetizations produced by several hundreds 
of electronic magnetic moments has recently been re- 
ported [4]. Attempting to reach the single spin resolu- 
tion, MRFM faces a number of experimental challenges. 
The observed reduction of the cantilever's quality factors 
for small separations between the tip and the sample [7,8] 
and fast spin relaxation [4] pose a problem in resolving 
a single spin signal. In this paper we analyze the latter 
effect. In particular, we study relaxation of a spin in the 
resonance slice [4] due to the thermal vibrations of the 
cantilever. We evaluate the spin relaxation rate and find 
a reasonable agreement with rates observed experimen- 
tally. We also propose that by shape engineering of the 
cantilever it is possible to filter the high frequency noise 
thus reducing the spin relaxation. 
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FIG. 1. Model setup for MRFM. 



In this work we consider the experimental setup used 
in Ref. [4]. One end of the cantilever is fixed, while the 



other end, to which micron size magnetic particle is at- 
tached, is driven to oscillate in a direction normal to the 
surface of the sample; see Fig. 1. The sample is doped 
with spin 1/2 magnetic impurities. The smallest flexural- 
mode eigenfrequency of the cantilever will be denoted by 
wo and the amplitude of the oscillations of the tip by 
zq. In typical MRFM experiments luq is several kHz and 
zo ~ 0.1-30 nm. The motion of the cantilever's mag- 
netic tip produces an oscillating magnetic field B t i P (r, t) 
at the spatial position r at time t (the origin is taken 
as the equilibrium position of the tip) . There is an ad- 
ditional constant magnetic field Bo applied normally to 
the sample, and a linearly polarized microwave magnetic 
field Bi cos (<x> r fi) in a direction perpendicular to Bo- 

The magnetic resonance condition is achieved as fol- 
lows. Let us assume for simplicity that a single impurity 
spin is positioned at point r immediately below the can- 
tilever tip, as shown in Fig. 1. Then, the magnitude of the 
field produced by the magnetic tip at the spin's position is 



fl° ip (r) + V z B^(r)z(t), where B° p (r) is the 



BtipM) 

field of the tip in equilibrium, z(t) is cantilever's displace- 
ment, |r| ^> z(t). Let the frequency w r f of the microwave 
field be chosen such that gfiB\B® ip (r) — Bo\ ~ hu) T i, where 
g is electronic g- factor and [1b is Bohr magneton. Then, 
in the rotating reference frame that rotates with fre- 
quency lu t { around the normal to the surface [9] , the spin 
experiences an effective magnetic field (see the inset to 
Fig. 1) 



B off = V,B& p (r)z(t) 



Bi/2. 



(1) 



As a result, the spin precesses around this effective field, 
thus following the direction of B or j, as well as cantilever's 
displacement, provided the latter varies with time adia- 
batically, i.e., luq <C w er r = g/is|B c ff|/?i. Thereby, the 
spin direction follows the cantilever motion and the spin 
exerts an oscillating force on the cantilever's tip with the 
amplitude, F = fiB^zB^ ip . In Ref. [4] the shift of the 
cantilever's frequency Awo due to the magnetization of 
the resonant spins has been measured. The signal that 
corresponds to roughly 100 fully polarized spins was ob- 
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served to decay on a time scale of 100 ms, thus indicating 
that the induced magnetization of the resonant spins re- 
laxes due to magnetic fluctuations or "noise" whose ori- 
gin we will discuss. 

In this letter we show that the main source of the noise 
causing spin relaxation is likely to be related to the can- 
tilever's thermal vibration. The cantilever displacement 
can be written as z(t) = z c (t)+Sz(t), where the first term, 
z c (t) = z cos (uj t) +zi, is due to the regular (driven) os- 
cillation of the cantilever tip, while 5z(t) is due to the 
thermal motion of the tip. Here, z\ defines the relative 
position of a spin within the resonant slice {\z\\ < z ). 
Substituting the above z(t) into Eq. (1) we obtain an 
effective Hamiltonian for a spin in the rotating frame 



H = H (t)+n{t)s z 



(2) 



where H = gp B {\7 z B° ip z c (t) s z + (B 1 /2)s x ], s is elec- 
tronic spin, and n(t) = #/z_bV z B^ ip Sz(t). In the follow- 
ing analysis we set g = 2 and the Planck constant % = 1, 
unless stated otherwise. We assume that the noise n(t) 
is Gaussian with the correlation function K(t\ — t 2 ) — 
{n{t\)n{t 2 )) to be specified below. We derive the Bloch 
equations for the spin from the dynamics governed by the 
Hamiltonian (2). We introduce the Keldysh contour and 
define a real time partition function along the contour as 



Z = T C exp 



Hdt exp 



Hdt 



(3) 



where T c denotes ordering along the contour [10]. Thus 
points on the forward branch (— oo — > oo) are ordered 
with increasing times, while points on the return branch 
(oo — ► — oo) are ordered with decreasing times. The su- 
perscripts / and r on i will indicate to which contour t be- 
longs. In Eq. (3) the time ordering operator T c sets oper- 
ators on the return branch in front of the operators on the 
forward branch of the contour. The expectation value of 
an operator 0{t) (O — s x ,s y ,s z ) can be obtained by us- 
ing the partition function Z as (0(f)) = Tr \T c O(tf)Z], 
where the trace is taken over the states of the spin and 
over the distribution of the classical noise n(t). Aver- 
aging the partition function Z over the noise under the 
assumption that the noise is Gaussian, we obtain 



/cr \ <T -i f °° H dt -if"" H dt 



x exp 



1 f°° 

-- K{t 1 -t 2 )s a z {t2)s a z {t 1 )dt 1 dtz 

* J -co 



, (4) 



where s z (t) — s{{t) — s r z (t). The effective action for the 
spin defined by Eq. (4) is generally non-local. It can be 
significantly simplified by introducing the Bloch-Rcdficld 
approximation [9]. If the interaction of the system with 
the noise is weak, one can expect that on a scale of the 
noise correlation time, r c , the spin's evolution will be 
mostly determined by the unperturbed Hamiltonian H - 



That is, in the non-local part of the action one can replace 
s a z (t 2 ) by exp[i H (T)dT]s a z (h)exp[-i £ H ( T )dT}. 
The latter operator can be evaluated assuming that 
Hq(t) varies adiabatically. By substituting the resulting 
expression into Eq. (4), the non-local part of the effective 
action in Eq. (4) can thus be approximated as 

AS = - i J | cos 9 [cos 9s a z (t) + sin 9s a x (t)} S n [0] 

+ sin 9 [sin 9s a z {t) - cos 9 s a x (t)]S n [u cS {t)]}dt , (5) 

where S n [uj] — K(t) exp (iu>t)dt is the power spec- 
trum of the noise n(t), to e s(t) = gpB{[^ zB^ ip z c (t)] 2 + 
_B^/4} 1 / 2 , and sin# = g^ B Bi/2uj cS {t). In the derivation 
of Eq. (5) we have assumed that w e ff / (duj e s / dt) 3> r c . 
The effective action in Eq. (5) is now local and allows for 
a straightforward derivation of equations of motion for 
the spin components s XtVtZ . Using s XtVtZ = i[H c s , 8 x , y , z ], 
where H c s is defined by Eqs. (4) and (5), and ordering 
the operators according to Eq. (3), we obtain a set of 
Bloch equations for the magnetization components 

s x = -w e s(t) cos 9(t)s y + f3(t)s z - a(t)s x , (6a) 
s y = uj cS (t) cos 9 (t)s x - w e s(t) sinO(t)s z - a(t)s y , (6b) 
s z = Lo cS {t) sin9(t)s y . (6c) 

In Eqs. (6), terms with coefficients a and (3 de- 
scribe relaxation of the magnetization due to noise, 
a(t) = {cos 2 0S„[O] + sm 2 9S n [uj cS (t)}}/2, and (3(t) = 
sin26'{S'„[0] - S , „[w e ff(t)]}/4. By solving the Bloch equa- 
tions in the adiabatic limit, i.e., by assuming that the 
coefficients in Eqs. (6) are slowly varying functions of t, 
we find that the component of the magnetization parallel 
to the effective field decays as 

|«eff(t)| - exp { - 1 J sin 2 9(T)S n [uj cS (T)}d T } . (7) 

Now we turn to the evaluation of the power 
spectrum of the magnetic noise produced by the 
cantilever tip. The magnetic noise is related 
to the mechanical noise of the tip as S n [ui] = 
(5MsV z B t ° ip ) 2 dtexp(iurt)(5z(0)6z(t)). It should be 
emphasized that we are interested in the mechanical noise 
of the cantilever at frequencies very high (of order u e ff) 
compared to the lowest cantilever eigenfrequency u> . As 
a consequence of the linearity of mechanical oscillator the 
high frequency modes of the cantilever will be essentially 
unaffected by the driving force at frequency ujq] rather it 
is the thermal driving force that is important. 

The energy of the vibrating cantilever can be expressed 
as E c = f dx[p{d t zf + EI(d 2 x z) 2 ]/2, where z(x,t) is 
the transverse displacement of the cantilever at point 
x and time t, p is cantilever's linear mass-density, E 
and I are Young modulus and moment of inertia of the 
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cantilever's cross-section respectively, and I is the can- 
tilever's length. The equation of motion for the free can- 
tilever is pd 2 z = Eld^z [11]. The general solution to 
this equation can be written as z(x, t) — ^ n z n (t)(j) n (x), 
with tfin's satisfying the eigenvalue equations for u> n : 
pu) 2 4> n — Eld^n, supplemented by appropriate bound- 
ary conditions [11]. In the eigenfunction basis 4> n , the 
energy of the cantilever can be expressed in terms of the 
amplitudes z n : E c = J2 n [ rn ^n+ rnul n z n\/^j where m = pi 
is cantilever mass. Hence, the cantilever can be mod- 
eled as a collection of independent oscillators. In ther- 
mal equilibrium we then have that (5z„ (ui)Sz n ' (w)} = 
S n ,n' (^ksT / mu>^)[6(u) n + u>) + 5(u n — u)], where fc^T is 
the temperature and the delta functions S(u> n ± u>) have 
effective width of order of damping coefficient for the n- 
th mode. The magnetic noise power spectrum can be 
expressed as 

SM = ( 3MB V z 5 t °ip) 2 E ^(lWznH\ 2 ) , (8) 



where the 



cigenfunctions (f> n (x) are normalized, 
jQ(j)^ l {x)dx = I. The power spectrum given by Eq. (8) 
generally represents a set of distinct peaks. However the 
effective frequency of a spin, w e ff(i), sweeps over a large 
number of peaks. Moreover, the form of Eq. (7) indicates 
that it is the integral characteristics of the spectrum that 
are relevant for the determination of the spin relaxation 
rate. Therefore the sum in Eq. (8) can be replaced by the 
integral, ^2 n — ► J (dn/ 'duj n )duj n and the finite width of 
the delta functions <5(u>„±u;) can be neglected. Then, tak- 
ing into account that for n ^> 1, u) n = (nn/l) 2 (EI/ p) 1 ^ 2 , 
4>l(l) = 4, and uj q = (3.52/l 2 ){EI / p) 1 ' 2 , sec Rcf. [11], 
Eq. (8) yields the average noise spectral density 



(g,i B V z B°) 2 k B T 



mujo 



1/2 5/2 



(9) 



Substitution of Eq. (9) into Eq. (7) gives the spin re- 
laxation rate. For large times, t ^> uj^ 1 and for 
\7 z B® ip \z ± Zi\ Bi, the integral in Eq. (7) can be 
readily evaluated yielding |s e ff(i)| ~ exp (— t/r m ), where 
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(10) 



Here we have introduced spring constant k s = muj 2 and 
LUl = p B Bi/h. 

We can now compare the relaxation time given by 
Eq. (10) with the experimental results of Ref. [4]. For 
tip-sample separations of 800 nm the gradient of the mag- 
netic field (V z B® ip ) was found to be roughly 1 x 10 5 T/m. 
Taking other experimental parameters from Ref. [4], we 
obtain from Eq. (10) the maximum T m ~ 350 ms (for 
z\ =0). Spins located away from the center of the res- 
onating region will have shorter relaxation times. This 



as well as the simplifying theoretical assumptions about 
the cantilever geometry are likely to be responsible for 
about five times difference between the theoretical esti- 
mate for the longest relaxation time and the experimental 
results [4]. We also note that because of the relaxation 
time distribution, measurement of an ensemble of spins 
is likely to yield a non-exponential MRFM signal decay. 

In the mechanism discussed above, the spin relaxation 
is caused by the thermo-mechanical noise of the high fre- 
quency modes of the cantilever tip. This suggests a way 
to reduce the relaxation rate by engineering the shape 
of the cantilever to reduce this noise. The simplest ap- 
proach that we analyze here is to place a massive particle 
(M) on the tip. This will effectively filter out frequencies 
uo n for which n > m/M. For such high-frequency modes 
it is easy to show that 4> 2 n {l) ~ 2yJ p 3 EI/(M 2 cj n ). For 
example, for a particle with M — m/10 and the same 
experimental parameters, the relaxation time would in- 
crease to t ~ 5 s. 

In summary, we have found that high frequency 
thermo-mechanical noise associated with high frequency 
modes of the cantilever can induce significant spin relax- 
ation in the magnetic resonance force microscopy, and is 
a likely reason for the coherent signal loss in the recent 
high-sensitivity experiments [4]. To reduce the influence 
of this noise, we have proposed cantilever shape engineer- 
ing, which can lead to significant enhancement of MRFM 
sensitivity. 
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